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The Quantum Spin Hall insulator is characterized by the presence of gapless helical edge states 
where the spin of the charge carriers is locked to their direction of motion. In order to probe the 
properties of the edge modes, we propose a design of a tunable quantum impurity realized by a 
local gate under an external magnetic field. Using the integrability of the impurity model, the 
conductance is computed for arbitrary interactions, temperatures and voltages, including the effect 
of Fermi liquid leads. The result can be used to infer the strength of interactions from transport 
experiments. 



The quantum spin Hall effect (QSHE) is a property 
of certain two-dimensional electron systems with strong 
spin-orbit coupling [TJ [5] . The bulk of the system is elec- 
trically insulating, while a conducting "helical edge" ex- 
ists at the boundary in which electrons of opposite spin 
move in opposite directions [3- 5 . Due to this reduction 
of the number of degrees of freedom, the QSHE edge is 
expected to realize the physics of a spinless Luttinger liq- 
uid, as opposed to a conventional one-dimensional wire 
that represents a spinful Luttinger liquid [BJ. The Lut- 
tinger liquid is the generic state of metallic interacting 
electrons in one dimension [7] , while metallic electrons in 
higher dimensions typically form a Fermi liquid. 

The QSHE is realized in (Hg,Cd)Te quantum wells 
[9] where measurements of the conductance indicate the 
existence of helical edge modes. The simplest measure- 
ment to perform on such a system would be a two- 
terminal conductance measurement. Such a measure- 
ment can confirm that the current is carried by helical 
one-dimensional edge channels, but it can neither pro- 
vide information on the interaction strength within those 
channels, nor verify the expected Luttinger liquid behav- 
ior. This is the case because, when a clean interacting 
wire is placed between Fermi liquid contacts (modeled 
as a non-interacting wires) , the measured conductance is 
insensitive to the interactions [TOl [H] . 

As it turns out, there is a way in which the two ter- 
minal conductance can provide information on the in- 
teraction strength within the edge modes. A common 
way of studying one-dimensional systems, both theoreti- 
cally and experimentally, is by exploring impurity effects 
on measurable quantities such as their conductance. In 
general, the problem becomes quite involved when in- 
teractions are present, and one usually has to rely on 
the asymptotic behavior of such quantities (at high or 
low temperatures, for example) to extract information 
on the interaction strength. However, in some unique 
cases certain properties of the edge model make it pos- 
sible to obtain exact solutions. The QSHE edge is an 



example of such a system, since the model of a spinless 
Luttinger liquid with an impurity is "integrable" |12j . 

In order to utilize the powerful tool of integrability to 
describe actual measurements on a QSHE edge, backscat- 
tering must be induced within a single edge (the model 
describing backscattering between the two edges of the 
QSHE system is not integrable). In principle, this can be 
done by means of a magnetic impurity that locally breaks 
time-reversal symmetry. However, it is much more desir- 
able to find a way to engineer an impurity with a tunable 
strength, in order to induce the crossover between weak 
and strong backscattering. 

In this work we consider combining the effects of an 
externally applied magnetic field and a local gate voltage 
to form an artificial impurity on the QSHE edge. The 
magnetic field direction is carefully chosen such that it 
breaks time-reversal symmetry yet leaves the edge modes 
gapless. These edge modes, now unprotected, become 
sensitive to the local perturbation generated by the gate 
in the form of an induced Rashba spin-orbit coupling. 
The strength of the impurity is set both by the mag- 
netic field and the gate voltage. With controlled means 
for introducing an impurity, the integrality of the edge 
model [T21fbl] allows us to extract the shape of the non- 
equilibrium, finite temperature conductance curve, which 
strongly depends on the value of the Luttinger parame- 
ter. Hence measuring the conductance throughout the 
crossover from weak to strong backscattering could pro- 
vide information on the interaction strength within the 
edge channels. 

The setup we have in mind (see Fig. [TJ is similar in 
spirit to a quantum point contact in fractional quantum 
Hall effect (FQHE) devices [TS]. There, particle backscat- 
tering between modes with opposite chirality is enhanced 
with the aid of two gates depleting the electron density 
and bringing the two edges of the sample closer together. 
However, for the QSHE device we consider, backscatter- 
ing between counter-propagating modes takes place on 
the same edge. Hence, we do not require that the two 
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FIG. 1: Schematic diagram of the proposed experimental 
setup. Voltage the top gate is used to locally tune the strength 
of the Rashba spin-orbit coupling. Combined with a magnetic 
field aligned along the electron spin quantization axis, a gap 
appears in the edge spectrum, giving rise to backscattering in 
the helical edge. 



edges of the sample be brought together, and a single 
gate is sufficient. Recently, leading corrections to the lin- 
ear conductance induced by a generic magnetic impurity 
in a fractional topological insulator were calculated [IB] , 
There, unlike the integer case we study, an edge with 
repulsive interaction can be stable to magnetic perturba- 
tions. 

Note that although both the QSHE and the FQHE 
edges realize a spinless Luttinger liquid, the Luttinger pa- 
rameter for the QSHE can in principle obtain any value, 
while for the FQHE it is restricted to quantized values. 
Another crucial difference between the two systems is 
embodied in the effect of Fermi liquid contacts discussed 
earlier. For the FQHE, contacts are expected to have no 
effect on the conductance, due to the spatial separation 
of modes of opposite chirality. This has been observed 
in experiments [T7HT5] . Therefore, the QSHE case has 
the potential to provide the first experimental test of in- 
tegrability at non-quantized values of the Luttinger pa- 
rameter and in the process verify the effects associated 
with Fermi-liquid contacts. 

We start by considering the non-interacting case, solv- 
ing the scattering problem of two gapless regions sepa- 
rated by a finite strip in which an energy gap is present. 
We find the reflection strength and show that it can dis- 
play resonant behavior for some values of the parameters. 
We then consider interactions and use a method known 
as the thermodynamic Bethe ansatz to obtain the non- 
equilibrium finite temperature conductance for various 
values of the Luttinger parameter [TJ] . 

The low energy physics of the non-interacting edge in 
the presence of a magnetic field B and a position depen- 
dent Rashba spin-orbit coupling a(x) is described by the 



Hamiltonian 



H = —ihvFa z d x 
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where vf is the Fermi velocity, the cr's are the Pauli ma- 
trices, {., .} denotes an anticommutator, g e is the electron 
Lande ^-factor and fig is the Bohr magneton. To sim- 
plify notation, in the following we take h — 1 and define 
M = fiB9eB/2. For M = a = the spectrum of this 
Hamiltonian is gapless, E = ±vfP- When a magnetic 
field is turned on, the energy spectrum becomes gapped, 
unless the magnetic field is parallel to the spin quantiza- 
tion axis of the electron. In that case the effect of the field 
is merely to shift the Dirac point and E = ±(vfP + M). 
In the absence of a magnetic field, a finite constant spin- 
orbit interaction a(x) = ao renormalizes the electron ve- 
locity to v a = \J 'oq +v F , and rotates the electron spin 
quantization axis by an angle cos 9 — vp/v a about the x 
axis [2D] . Note that the spins of the counter-propagating 
modes remain anti-parallel in the presence of the Rashba 
term as required by time-reversal symmetry. 

Let us now consider a system in which the magnetic 
field is uniform and points along the spin quantization 
axis, while a finite constant Rashba coupling exists only 
within a finite strip of width d, a(x) = ao<d(x)<d(d — x). 
Outside the strip (x < 0,x > d), the energy spectrum is 
gapless, while within the strip (0 < x < d), the external 
magnetic field is no longer aligned with the spin polar- 
ization axis, and the energy spectrum becomes gapped 



E = ±yJ (vlp + v F M/v a f + alNP/vl, 



(2) 



with the energy gap E g = \2a$M jv a \. In the presence 
of the external field, the two otherwise decoupled spinors 
now mix in the region combining both the field and the 
spin-orbit coupling. The result is a square scattering bar- 
rier, from which incoming waves can be reflected. In 
the limit of a narrow constriction, this region acts as a 
localized impurity in our helical quantum wire, whose 
strength is controlled by M and ao. In reality this can 
be realized by varying the voltage of a nearby electro- 
static gate which enhances the Rashba coupling in the 
vicinity of the gate, while the Rashba coupling far from 
the gate is negligible [21] . 

We solve the scattering problem by defining the scat- 
tering state in each region to be 



ip R e ipRX + rip L e lpLX x < 

a + ijj + e ip + x + a-tp-e ip - x < x < d 



> d 



where ipR — (L0), ipL — (0, 1) and ip± — (iap±,VFP± — 
M - E). The momenta p R/L = (±E - M)/v F corre- 
spond to the right (R) and left (L) movers outside the 
strip, while p± are the two momenta inside the strip, 
corresponding to the solutions of ([2| at a given energy. 
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The nontrivial part of the solution for r and t, the reflec- 
tion and transmission amplitudes is to find the correct 
matching condition for the wave function ^ at x = 0, d. 

For a general profile of a(x), the Schroedinger equation 
Hip = Ei/j [Eq. 0] can be solved formally as ^(xi) = 
T Xl ,x ^(xo) where the transfer matrix is written as 

T XuX0 =P x e J *° 4+- 2 I - 21 »J, (3) 

with P K representing the path ordering operator. For a 
step in a, a(x) — a>oO(x), we set Xq — —6 and X\ = 8 
and then take the limit of 6 — >• 0. The contribution of the 
terms including the magnetic field and the energy in the 
exponent will vanish, and we are left with the matching 
condition 

/ \ 1/2 

^(0+)=f^J e l0oa ^(O-), (4) 

with tan 2^o = &o/vf- 

Using the boundary condition Q at x — 0, and a sim- 
ilar one at x ~ d, we obtain the solutions for r, t. The 
analytical form of these solutions is lengthy, therefore 
we will not present it here but rather plot the reflection 
probability R = \r\ 2 in Fig. [2J The behavior of R as a 
function of the field at a fixed (nonzero) value of cko/vf 
can be described as follows: at zero field, the reflection 
is zero, while for large fields M ^> v a E/cto, the reflection 
is perfect, R — 1. In between, R can display two types 
of behaviors, depending on whether the waves inside the 
barrier are evanescent or propagating. For evanescent 
waves E 2 < ot^M 2 jv 2 a , R rises monotonically towards 
unity, while propagating waves result in an oscillating 
reflection amplitude, due to Fabry-Perot type of interfer- 
ence resonances. The condition for a resonance is simply 
(p+ — p~)d = 27m, and depends both on the value of 
ao and M. Such Fabry-Perot resonances have also been 
predicted for a QSHE edge state under a spatially inho- 
mogeneous magnetic field [22] . 

While the resonances appearing for the finite barrier 
would provide a test for the existence of helical edge 
modes in the absence of interactions, it is expected that 
interaction effects are important in ID systems and can 
renormalize drastically the backscattering created by a 
single impurity [53]. This is true even for weak repul- 
sive interactions since single electron backscattering is 
described by a relevant operator (in the renormaliza- 
tion group sense), leading to a crossover from a weak- 
backscattering regime to a strong backscattering regime 
as the temperature is lowered. The non-interacting solu- 
tion is still useful in estimating the bare backscattering 
strength and its dependence on ao and M. In our dis- 
cussion of the non-interacting problem we considered a 
region of finite width d as the scatterer. In order to cal- 
culate conductance in the interacting case we need to 
consider a point like scatterer. The bare backscattering 
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FIG. 2: Reflection probability R as a function of the nor- 
malized magnetic field M for various Fermi energies E, and 
oto/vF = 0.1. The energy unit is Eo = hvF/d. 



strength for such an impurity could be estimated from 
our previous calculation by taking the limit of a very 
narrow barrier |30j . For weak magnetic fields and small 
ao it is simply given by R ~ (M ot2k / v F ) 2 , where a.2k is 
the 2/c component of the Fourier decomposition of a(x). 

The Hamiltonian of the interacting QSHE edge in the 
bosonization language is 

H= i^ g j d <9M 2 + (d x <t> L ) 2 , (5) 

where 4>r/l are left moving and right moving boson fields, 
g is the Luttinger liquid parameter and v is the edge ve- 
locity renormalized by interactions [24] . A backscattering 
term couples to 4>l — 4>r 

flB = Acos(0i(O)-0«(O)). (6) 

By defining even and odd non-local combinations of the 
fields C / O = 1/v2(0l(x, t)±cj)fi(— x,t)), the backscatter- 
ing term couples only to <j> , and the Hamiltonian breaks 
into two decoupled contributions. The part of the Hamil- 
tonian describing the odd fields is integrable, since it 
is identical to the massless limit of the boundary sine- 
Gordon (SG) model \12\ EH]- The even field theory is 
free and does not interact with the impurity. 

The integrability of the SG model was previously used 
by Fendley, Ludwig, and Saleur (FLS) OH [T3] to calcu- 
late the non-linear conductance in a point contact geom- 
etry for fractional quantum Hall states at v = 1/m. In 
Ref. [25] a similar formula is derived for the conductance 
at v = 1 — 1/m by exploiting the relation between the 
Kondo problem in a magnetic field and the SG model at 
finite voltage. Here, we use the same method to compute 
the differential conductance curves at m = 3, 4, 5 in order 
to demonstrate the behavior of the QSHE edge transport 
atl/2<<7<l. A result for continuously variable g can 
in principle be computed using a more involved technique 
developed in (5^]. Also, if the repulsive interactions are 
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weak (g close to 1), a more direct approach based on re- 
summed perturbation theory [27] can be used, but it is 
currently unclear whether interactions are weak in the 
QSHE edge. 

For g = 1 — 1/m, the current along the QSHE edge is 
given by [2"5] : 



I(V,T B ,T) 



In 



T(m-i) r de x 

2 7 cosh 2 [(9- ln(T B /T)] 

1 + e (m-l)V/2T-e+(0) 1 + g-(m-l)V/2T-6 + (oo) 
1 + e -(m-l)V/2T-e + (0) 1 + e (m-l)V/2T- e+ (oo) 



(7) 



Here Tb is an energy scale related to the impurity 
strength A by T B = CA 1 ^ 1 "^ [HES], where C is a non- 
universal (cutoff dependent) constant, e+ is the quasi- 
energy of the kink solution of the sine-Gordon model, 
and 9 the rapidity. The energy e+(9) of the kinks is 
computed numerically by solving a set of coupled inte- 
gral equations obtained from the thermodynamic Bethe 
ansatz. The full details of these equations are given in 
supplementary material. 

To account for the effect of non-interacting leads on 
the calculation of the conductance, we adapt a result 
from Ref. 123 where a self-consistency condition was de- 
rived for the chemical potential of the various excitations 
inside the wire, which is not equal to the external applied 
voltage. The consequences of this self-consistency condi- 
tion were extensively explored in Ref. [55J for g = 1/m. 
Here we carry out a similar analysis for g = 1 — 1/m. 
Denoting the chemical potential for the kinks and anti- 
kinks by n± = ±(to — l)W/2 and the external voltage by 
V, the self-consistency condition for W is [25J 



v = ~ (l- - ) /(H i • ir. 



(8) 



The results for the differential conductance G = dl/dV, 
with and without the contact correction, are presented 
in Fig. [3] The asymptotic behavior of G as a function 
of Tb/T matches the known predictions [23] . namely 
G ~ e 2 /h(T/T B ) 2/9 ~ 2 at low temperature, and G ~ 
e 2 /h (l - (Tb/T) 2(1 ~ 9) ") at high temperature. 

Though we have computed the full curve for particular 
values of g, all curves show similar features and a compar- 
ison with experimental data should confirm the expected 
Luttinger liquid behavior and yield a good estimate for 
the value of g. Note that when contact corrections are in- 
cluded the conductance always saturates to e 2 /h at high 
voltage or in the absence of backscattering. Nevertheless, 
the curve shape itself highly depends on g, and in par- 
ticular, the exponents of the asymptotic behavior remain 
the same as without the correction. 

The feasibility of our proposal depends both on the 
stability of the QSHE edge in presence of a magnetic 
field and the spin direction of the modes. The behavior 
of the QSHE under a magnetic field has been studied in 
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FIG. 3: Differential conductance for different values of g = 
1 — 1/m. Top figure: G = dl/dV (with contact corrections) 
as a function of Tb/T. Bottom figures: G = dl/dV as a 
function of V/T for two different values of T B /T. The differ- 
ential conductance has been scaled with the value Goo = e 2 /h 
and (1 — l/m)e 2 jh with and without contact corrections, re- 
spectively. The curves without contact corrections have been 
shifted down by 0.3 for clarity. 



an experiment where the conductance was measured for 
various tilt angles of the field with respect to the plane of 
the 2D electron gas [H [5J [2H] • The results show that on 
top of the contribution from the Zeeman coupling, when 
the field is perpendicular to the plane, the conductance 
drops rapidly with the field strength due to orbital effects. 
Nevertheless, a peak in the conductance of typical width 
B = lOmT exists at T = 30mK. Orbital effects result in 
an effective g-factor values of 20 — 50, the typical Fermi 
velocity is estimated to be vf = 5.5 x 10 5 m/s, and ao ~ 
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5 x 10 4 to/s. Therefore, even under the most restrictive 
conditions one can obtain a gap size of E g 100— 300mK 
in the vicinity of the gate in our setup. 

In our analysis we have made the assumption that the 
spin quantization axis far from the gated region is fixed 
along the edge, and therefore it is possible to align the 
magnetic field such that it does not gap out the edge 
modes in those regions. In principle, the preferred spin 
quantization axis along the edge is determined by the 
properties of the material, is not protected, and may tilt 
in a complicated way along the edge due to fluctuations 
of the Rashba coupling. However, if the edge is made 
smooth enough it is reasonable to assume that fluctua- 
tions in the Rashba coupling have a much smaller effect 
on the edge states than the intentional coupling induced 
by the gate, and therefore our analysis remains valid. 

To summarize, in this paper we proposed a realistic 
and controlled setup for transport experiments on the 
QSHE edge using the combined effect of time-reversal 
symmetry breaking and induced spin-orbit coupling. Us- 
ing the integrability of the resulting edge model we pre- 
dict the form of the non-linear conductance, which could 
be compared with future experimental data. 
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Supplementary material: Thermodynamic Bethe 
Ansatz equations 

The thermodynamic Bethe ansatz equations relevant 
for our particular model are the following coupled inte- 
gral equations [25] : 

L k (6') = In (l + e -(^( 9 ')+M*)/^ , (A.9) 

In this set of coupled integral equations, the indices j and 
k label the different excitation of the SG model. There 
are three types of excitations: breathers (1, ..,m — 2), 
kinks (+) and the antikinks (— ). The functions e k label 
their energies (parametrized by an angle 9 known as the 
rapidity). Kinks and anti- kinks are the charge carrying 



excitations of positive and negative charge respectively, 
while the breathers are neutral bound states of a kink 
and an anti-kink, and therefore their chemical potential 
is alway set to zero. In general, e + = e_, and for v ~ 
1 — I/to, the matrix element Nj k = 1 if the nodes j and 
k are connected in the following diagram 




otherwise Njk = 0. 

We solve these equation numerically for arbitrary val- 
ues of m. Then the energy of the kinks is used in Eq. ^ 
to compute the conductance without contact correction, 
or combined with the self-consistency condition defined 
in Eq.Q, to evaluate the conductance in presence of con- 
tact corrections. 



